Exact solution of a class of critical dynamical systems: Information routing in 

complete graphs 
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It is a controversial question to which extend scaling is observed in critical dynamical systems. 
Here we examine a class of dynamical systems describing the propagation of conserved quantities, 
such as the routing of information packages, and which are designed as critical systems. The long 
term information flow is governed by properties of the cyclic attractors in phase space. 

We consider two classes of information routing models, Markovian routing without memory and 
vertex routing involving an one-step routing memory. Investigating the respective cycle length 
distributions for complete graphs we present analytic approximations which become exact in the 
limit of large numbers of vertices, the thermodynamic limit. We flnd log corrections to power-law 
scaling for the mean cycle length, as a function of the number of vertices, and a subpolynomial 
growth for the overall number of cycles. 



Introduction." The propagation of perturbations is a 
central notion in dynamical system theory. One speaks 
of a frozen state when a perturbation tends to die out, on 
the average, during the course of time evolution and of a 
chaotic state when perturbations tend to spread out [l|, 
|2| . A given class of dynamical systems may change from 
frozen to chaotic behavior as a function of parameters, 
being critical right at the transition point. 

At criticality, information is on the average conserved, 
as one can regard a perturbation of a state as the in- 
formation about the persistence of small differences. A 
well studied example of a critical dynamical system is the 
Kauffman net with connectivity K ~ 2, an example of a 
random boolean network (3|, |j] . In statistical mechanics 
critical systems are generically scale invariant [5|, and it 
has been widely assumed that this statement would also 
hold for critical dynamical systems. Indeed numerical 
simulations seemed to support scaling in critical boolean 
networks, notably a vN scaling for the number of at- 
tractors as a function of the number of vertices TV had 
been proposed [3|, |j] . 

An important clarification then came with the exact 
proof that the number of attractors actually grows faster 
than any power of N , and that the results of the numer- 
ical simulations suffered from systematic undersampling 
of phase space [6|. It could be shown, on the other side, 
that the number of frozen and the number of relevant 
nodes in a large class of critical boolean networks obeys 
powerlaw scaling [7| . 

The situation is then that certain properties of criti- 
cal dynamical systems, at least for the case of random 
boolean networks, obey powerlaw scaling while others 
do not. It is hence important to investigate the possi- 
ble occurrence of scaling in different classes of dynamical 
systems. We study here a novel class of dynamical sys- 
tems modelling the transport of conserved quantities in 
network structures and which are hence critical per con- 
struction. 

Transport on networks, like the spreading of rumors [8| 



and diseases [9| in social networks or the flow of capital in 
financial networks [lO| is a basic process in biology 11| . 
as well as in sociology and technical applications. 

We considered here the transport of a conserved quan- 
tity, like information packages, in terms of routing pro- 
cesses. Information packages are sent from node to node 
and are routed at every vertex, as illustrated in Fig.[T] All 
routing process eventually ends up in one of the cyclic at- 
tractors. It has been shown previously that the geometric 
arrangement of the attractors on the network gives rise 
in the thermodynamic limit to a non-trivial distribution 
for the information centrality, which measures the num- 
ber of attractors intersecting at a given vertex [ij] . We 
present here analytic solutions for two types of vertex 
routing models, Markovian routing in the absence of a 
routing memory and vertex routing in the presence of an 
one-step memory. 

The analytic solutions are asymptotically exact in the 
thermodynamic limit A^ — > oo, they can be evaluated for 
large networks containing thousands to millions of sites. 
We present results for the scaling behavior of the overall 
number of attractors and for the mean of the cycle length 
distribution. 

Models.- The two classes of models we consider differ 
with respect to the absence/presence of a routing mem- 
ory. The phase space volume n is respectively linear and 
quadratic in the number of vertices N. 



• For the Markovian model the selection of the 
next active vertex is independent of the previous 
state [13|. At every point in time only one vertex 
is active, the vertex with the information package. 
The phase space is hence identical with the collec- 
tion of vertices; il ^ N; 

• For the vertex routing model the phase space is 
given by the collection of directed links; il = 
N{N — 1). At every point in time one directed 
link is active, the link currently transporting the 
information package, compare Fig. [1] 
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FIG. 1: Vertex routing dynamics for a A*' = 4 complete graph 
(a) A realization of the routing tables. Routing through the 
first vertex follows T312 = T213 = 1214 = 1, with all other Tnj 
vanishing. There are three cyclic attractors, namely (123), 
(243) and (1342). (b) Enumeration of all N(N - 1) = 12 di- 
rected edges, the phase-space elements, (c) The phase-space 
graph, (d) The same realization of the routing table as in 
(a), now in terms of the phase-space graph. 



In both setups the routing of information packages is 
realized through static routing tables. For every incom- 
ing edge the routing table specifies an allowed outgoing 
edge. A vertex k will transmit an information package, 
which was received from a vertex j, to a specific neigh- 
boring vertex i. The vertex routing table T corresponds 
to a tensor of binary elements Ti^j = {T)ikj G {0, 1}, 



T, 



ikj 



no routing allowed 

1 routing from Cji,. to Cki 



(1) 



where ej^ denotes a directed edge from vertex j to vertex 
k. An example of a routing tabic for a four-site network 
is presented in Fig. [1] 

We consider here critical models, viz models where the 
number of information packages is conserved. When the 
information is received along edge Cjk, it can hence be 
transmitted along only one outgoing edge Cki, 



^Tikj = l, ^T^kj - 



Zk 



(2) 



which is drawn randomly. Here z^ is the degree of vertex 
k, which is TV— 1 for fully connected networks considered 
here. For the Markovian model the routing table Ti^j is 
independent of j . 

The two different routing models represent two limit- 
ing cases for the modelling of rumors in social networks. 
Markovian dynamics describes the spreading of rumors 
when the spreading is independent of social relations be- 
tween the sender and the receiver. Vertex routing pro- 
cesses are, on the other side, important when the origin of 
the rumor is taken into consideration. A person receiving 
a work-related rumor might be prone to pass the rumor 
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FIG. 2: Random walks through configuration space for the 
Markovian model (left) and for the vertex routing model 
(right). In order to find an attractor independent of the size 
of their basins of attraction (light color) one needs to close 
the path at the respective starting points. 



to another colleague, rather than to a family member or 
a friend. 

Solution of the vertex routing model- The dynamics 
consists of random walks through configuration space, as 
illustrated in Fig. [2] One can hence adapt the consid- 
erations [4I, used for solving the Kauffman network for 
large connectivities i^ — ?► 00, in order to solve the vertex 
routing model analytically. 

At time-step i = 0, 1, ... the number of phase-space 
elements visited previously is t. The probability that the 
next site is a site previously visited is then t/{N — 1). 
For the trajectory to close the routing has to retrace the 
existing path, which happens with probability 1/ {N — 1), 
compare Fig. [51 

The relative probability pt of closing the path at time 
t is hence 



t 



1 



Pt 



N - I N -I 



(3) 



We define now with qt the probability that the path is 
not closed at timestep t. The probability of having an 
open path at t + 1 is then 



Qt+i = qt{l-pt) = ^iH 



{N - 1)2 -i 



i\ (N^iy 



with gi = 1, or 



qt+i = 



1 



{{N-ir-iy. 

{N-iyt ((7V-l)2-t-l)! 



(4) 



The average number {Cl)v{N) of cycles of length L e 

[2, (iV - 1)2 + 1] is then given by 



{Cl)4N) 



1 N{N~1) 

L (iV-l)2 '^'^^-^ 

1 N{N-1) 



(5) 



((^-1)^-1)! 

L {N- 1)2^-2 ((iv - 1)2 - i -h 1)! 



The factor 1/L takes care of overcounting, the factor 
N{N — 1) is the number of possible starting phase space 
elements, the factor 1/{N — 1)2 the probability to close 
the path exactly at the starting phase space element and 
qt=L~i is the probability that a path containing L = t+1 
directed links is open. Note that the path has to close 
at the starting phase space element and not at any other 
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FIG. 3: The cycle length distributions Pv{L), rescaled by 
log(f2), for the vertex routing model. The dashed line, 2/L, 
represents the large-A'^ and small-L limiting behavior. Inset; 
Log-Linear plot of the respective numbers of cycles, as a func- 
tion of the phase space volume Q. 
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FIG. 4: The mean cycle lengths {I/)t,,5 for the vertex routing 
with quenched dynamics (blue circles) and the vertex routing 
with on the fly dynamics (green diamonds) respectively, as a 
function of the phase space volume f2; log-log plot. The re- 
spective dashed lines are fits using a-|-6%/r2/log(n)4-c/log(n) 
and a! -\- b'y/Ti, with a, a' , b, b' , c being fit parameters. 



previously visited directed link. The probability to find 
an attractor would otherwise depend on the size of its 
basin of attraction, compare Fig. [5J The relation ([5]) is 
an approximation of the real cycle length distribution as 
it doesn't take into account corrections for self intersect- 
ing paths. These corrections drop however as 1/N and 
can be neglected in the thermodynamic limit. We denote 
with 



PviL,N), 



J2pv{l,n) = i 



the normalized cycle length distribution for vertex rout- 
ing with memory, which is obtained from ([S]) by normal- 
izing {CL)viN). 

The cycle-length distribution for the Markovian model 
Pjn{L,N), can be derived in an analogous fashion, the 
phase space is here the collection of vertices. 12 1 

Intermodel scaling relation- The number of cycles of 
length L e [2,N] for the Markovian model is [l2J 



(CL),„(iV) 



1 



1 



iV! 



L {N - 1)^ {N - L)l 



(6) 



TABLE I: Scaling with the number of vertices A'", for the 
number of cycles and for the mean of the cycle length distri- 
bution, respectively for vertex routing (v) and the Markovian 
(m) model, and for the two types of dynamics, quenched and 
on the fiy. Only relative quantities can be evaluated for on 
the fiy dynamics. 







quenched 


on the fly 


(v) 


number 
mean 


log(iV) 

iV/log(iV) 


A 


(m) 


number 
mean 


log (A) 
\/A/log(A) 


^/A 



Substituting A^ by (A^ 
N the scaling relation 



1 in ([5]) we obtain for large 



{CL)viN) ^ {CL)rn{{N^lf + l) 



(7) 



between the number of cycles of the vertex routing and 
the Markovian model, {Cl)v and {CL)m- 

Results- Evaluating numerically the analytic expres- 
sions (O and ([ni for the number of cycles we find clear 
evidence for a logarithmic growth for the total number of 
attractors {n)y^„i = J2Li^L)v,m, as a function of phase 
space volume J7, see Fig. [3l The total number of cycles 
hence grows slower than any polynomial of the number 
of vertices A^, in contrast to critical Kauffman models, 
where it grows faster than any power of A^ [6| . 

The normalized cycle length distributions Pv,m{L) 
thus scale as l/log(f7), due to the normalization. The 
rescaled distributions \og{n) p^,rn{L) approach the ther- 
modynamic limit rapidly, compare Fig. [3l For small cy- 
cle lengths L the limiting functional form of the rescaled 
distributions is 2/L, while for large L — > L^ax it falls 
off as LinH^x"' iQ-L-max^ ^jjg limiting behavior of 

log(r2)pt,,m(i) is identical for both models, due to the 
intermodel scaling relation ([7]). 

The determination of the scaling behavior is somewhat 
more subtle for the mean cycle length, see Fig. 01 We find 
that the functional dependence on the phase space vol- 
ume is best reproduced by a + 6a/51/ log(fi) -I- c/ log(r2), 
where a, &, c are free parameters, which fits the data by 
about one order of magnitude better than a pure power 
law Ansatz a' + h' Vf . An overview of the obtained scal- 
ing relations is given in Table HI Note that in Figs. [3] and 
|4]we present only the data for the vertex routing model 
as it completely overlaps for large phase spaces f2, due to 
the scaling ([7]) , with the results for the Markovian model. 

On the fly dynamics- In addition to working with pre- 
determined vertex routing tables one can generate dy- 



namics 'on the fly' without exphcitly creating routing 
tables for ah vertices of the network. 

For this kind of dynamics a routing for a given ver- 
tex is selected only when the trajectory visits this ver- 
tex. A cyclic attractor is then found when one state 
of the phase space (edge or node) is visited more then 
once. Then, one can show that the probability to 
find an L-cycle is pv{L,N) = J2t=L'' 1t^ where we de- 
noted with Pjj{L,N) the weighted cycle length distri- 
butions for the vertex routing model and where qt is 
given by (j4]) . An analogous relation holds for the Marko- 
vian model. Generalizing the scaling relation ([7|) one 
finds Pv{L, N) = Pm{L, {N — 1)^ + 1) and consequently 
{L)y{N) = {L)rh{{N - 1)^ + 1), where p„ denotes the 
weighted cycle length distributions for the Markovian 
model. 

We find clear evidence for a scaling ^ VTj for the mean 
cycle lengths, compare Fig. IH and no indications for log 
corrections. Note that the overall number of cycles can- 
not be obtained when routing on the fly, only relative 
quantities can be evaluated. 

Discussion and Conclusions- For boolean networks 
the phase space volume fi is 2^ and hence grows expo- 
nentially with the number of vertices A^. The fact p that 
the number of attractors grows faster than any power of 
A^ could hence in principle be related to the exponential 
growth of the phase space volume. Whether these two 
scaling relations, the exponential growth of the phase 
space volume and the superpolynomial growth of the 
number of attractors, are indeed causally related is how- 
ever not known to date. 

Our results indicate that scaling in critical dynamical 
systems may generically be non-universal. For two solv- 
able models with polynomial phase space scaling, as a 
function of the number of vertices, we find a subpoly- 
nomial growth in the number of attractors and log cor- 
rections to powerlaw scaling for the mean cycle length. 
These results indicate that there is probably no generic 
relation between the scaling of the phase space volume 
and the number of cycles. 

We also note that other properties of critical dynam- 
ical systems, like the scaling of the number of frozen or 
relevant nodes for critical boolean networks [7[ , may show 
highly non-trivial behavior. For the case of vertex routing 
models one may define a measure of centrality, informa- 
tion centrality, determined by the number of attractors 
intersecting a given vertex, which scales to a non-trivial 
limiting distribution in the thermodynamic limit 



12|. 



Our results may also be seen in the context of the surge 
in interests in modelling [14l [l5| and in experimentally 



investigating [ly, ll7| the spontaneous neural dynamics 
of the brain. The observation of powerlaw scaling rela- 



is generically dependent on the observation modus. The 
unbiased statistics of a certain property, like the number 
of attractors or avalanches, may differ from a statistics 
obtained via stochastic sampling {pv,m{L) and Pv,m{L) in 
our case). The later will in general be dependent on the 
size of the respective basins of attraction of the dynam- 
ical process considered, viz of a cycle or an avalanche. 
For the case of the vertex routing model studied here we 
found logarithmic corrections to powelaw scaling for the 
unbiased statistic and pure powerlaw scaling for stochas- 
tic on the fly sampling. 



tions 18| have been interpreted as evidence of a critical 
self-organized neural state [19| . Our results indicate that 
there is no universal relation in dynamical systems theory 
between criticality and powerlaw scaling and that scaling 
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